On improving single photon sources via linear optics and photodetection 
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In practice, single photons are generated as a mixture of vacuum with a single photon with 
weights 1 — p and p, respectively; here we are concerned with increasing p by directing multiple 
copies of the single photon-vacuum mixture into a linear optic device and applying photodetection 
on some outputs to conditionally prepare single photon states with larger p. We prove that it is 
impossible, under certain conditions, to increase p via linear optics and conditional preparation 
based on photodetection, and we also establish a class of photodetection events for which p can be 
improved. In addition we prove that it is not possible to obtain perfect (p = 1) single photon states 
via this method from imperfect {p < 1) inputs. 

PACS numbers: 03.67.-a, 42.50.Dv 



Single photon sources are important, for example to 
achieve secure quantum key distribution Q and for lin- 
ear optic quantum computation Q , yet generating single 
photons remains challenging. The traditional method in- 
volves photodetection on one output mode from a non- 
degenerate parametric downconversion process to post- 
select a single photon in the correlated mode More 
recently alternative single-photon sources have been em- 
ployed, including molecules 0, quantum wells 0], color 
centers |^ , ions Q and quantum dots Q . In these exam- 
ples, the probability of more than one photon being pro- 
duced is much lower than that for a Poissonian process, 
but the vacuum contribution can be quite high. Under 
ideal conditions, with a single mode output stable over 
time, the output state can be described by the density 
matrix 

p, = (l-p)|0)(0|+p|l)(l|, (1) 

with p the probability of obtaining a single photon; p is 
sometimes referred to as the efficiency for producing sin- 
gle photons. Whereas the theoretical single-photon state 
corresponds to p = 1, the experimental single-photon 
state is a mixture of the single photon with the vacuum. 

Increasing the value of p (the efficiency of produc- 
ing single photons) is an increasingly important objec- 
tive because of requirements for quantum optics experi- 
ments, especially those concerned with quantum informa- 
tion processing. Much of this effort is directed to improv- 
ing sources, but here we pose the question as to whether 
a source of single photons with efficiency p can be pro- 
cessed via linear optics and photodetection to yield fewer 
photons but with higher p. More specifically, would it be 
possible to direct N copies of these single photon states 
with efficiency p into an A^-channel passive interferome- 
ter (an interferometer that involves only passive, linear 
optical elements) to yield an output single photon with 



efficiency p' > p for certain photodetection records on 
the other outputs? If the answer is yes, then interferom- 
etry and photodetection could be employed to improve 
the efficiency of single-photon sources. 

In fact the answer will be no, provided we consider 
detection results where all but one of the photons are 
detected. This is the most straightforward way of en- 
suring that the output state contains at most one pho- 
ton. This restriction is necessary in most applications of 
linear optics quantum information processing since two- 
photon contributions can have unwanted side effects, e.g. 
allowing for certain multiphoton attacks in quantum key 
distribution applications. If we allow other detection 
results, it is possible for low-efficiency (small p) single- 
photon states to yield, via linear optics and conditional 
preparation based on photodetection, an output with a 
larger probability for a single photon. However, these 
schemes also yield non-zero probabilities for photon num- 
bers above one. 

A passive interferometer is comprised of beam splitters, 
mirrors, and phase shifters. Each of these elements pre- 
serves total photon number from input to output under 
ideal conditions. No energy is required to operate these 
optical elements, hence the term passive. These are also 
known as linear optical elements. More generally polar- 
ization transforming elements can be included, but here 
we are concerned only with a scalar field treatment; in 
fact polarization effects could be included by doubling the 
number of channels and treating the two polarizations in 
a mode as two separate channels. Mathematically, a pas- 
sive interferometer transforms the amplitude operators 
of the incoming fields a via the matrix transformation 
at ^ A^at with A S U(iV) 0. 

In the general case we start with a supply of N mixed 
states of the form For additional generality we al- 
low the different inputs to have different probabilities for 
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FIG. 1: Schematic setup of the network. We assume A*" 
incoming modes prepared in the state ^ with different pi. 



finding a single photon, pi, and we denote the maximum 
of these probabihties by Pmax- The initial state may be 
expressed as 
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where Ps = YiiPi'i^ ~ Pi)^ the vector s = 
(si,-- - ,sn)'^, {si — 0,1), gives the photon numbers 
in the inputs. The interferometer transformation i-^ 
A-^d^ results in the output state 
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Without loss of generality, we take mode 1 to be the 
mode in which we wish to improve the photon statistics. 
We perform photodetections on the other N — 1 modes, 
and examine the final state in mode 1 conditioned on 
the results of these photodetections. The total number 
of photons detected is D, and the number of photons 
detected in mode j is rij (see Fig. QJ. It is easy to see 
that no better result can be obtained by performing pho- 
todetections on fewer than — 1 modes; this would be 
equivalent to averaging over the photocounts for some of 
the modes. We assume ideal photodetection in this anal- 
ysis in order to determine the best results possible using 
linear optics and photodetection. 

The conditional state in mode 1 after photodetection 
in modes 2 to TV is 
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(4) 



ni=0 



where ni is the photon number in the remaining output 
mode (rather than the number of photons detected at a 
photodetector) . Each coefficient c„j is given by 



where |n) is a tensor product of number states in each of 
the output modes. The normalization constant K equals 



K = 



■ N 

E 

."1=0 



(JV) 



n 



(6) 



In order to find the expectation value in Eq. 10), we first 
introduce some notation. Let $s = {i\si = 1}, Eg = 
J2i Si (so Ss is the number of elements in and let 
Ys be the set that consists of all vectors comprised of the 
elements of $s- In addition we use the notation K' = 
K/ Y\f=2 • ^^'^ ~ ■ Using this notation, ci is 

given by 
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where 
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In order to determine if there is an improvement in 
the probability of finding a single photon, we need to 
determine the value of ci. However, determining ci 
requires evaluating K, which requires evaluating the 
expectation values in Eq. (0) for all possible values of 
ni. Instead, we will consider the ratio ci/cq. There are 
three main advantages to considering this quantity: 

1. The common constant K' cancels, so this expression 
is easily evaluated analytically. 

2. If ci/co is not greater than R = Pinax/(1 - Pmax), 
then it is clear that ci < Pmax- Thus we can determine 
those cases where there is no improvement. 

3. For Pmax <C 1, Co ~ 1 and R ~ Pmax- Therefore the 
improvement in ci/cq over R is approximately the same 
as the improvement in the probability of a single photon 

over Pinax- 

Ideally, we would determine the interferometer and de- 
tection pattern such that ci/cq is maximized, but this 
does not appear to be possible analytically. However, we 
can place an upper limit on ci/cq in the following way. 
Let us express the summation for cq as 
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(9) 



s;S,=_D+l fc:sfc 



where = Si except for = 0, and = rij except n[ — 
0. The quotient oi N — D takes account of a redundancy 
in the sum. Each alternative input s'^ has N — D zeros, 
so there are N ~ D possible alternative s that give the 
same s'^. We may reduce this quotient slightly if we take 
(5) account of the possibility that some of the inputs have 
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zero photon probability. Let there be iV — M inputs with 
Pi = 0, so the maximum total number of photons is M. 
If we limit the first sum in Eq. @ to s such that Pg 0, 
then the redundancy is Al — D. Therefore we obtain 



and (S'fc) to be as close to parallel as possible. For this, 
we consider the interferometer given by 



A21 = -e, A22 = Vl 
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M-D ^ ^ Pk 

s;Ps^O k;Sk = l 

where Sk = Sgk „/. Since we have limited the sum to 
terms where Pg 7^ 0, pk is nonzero, and thus the ratio 
(1 — Pk)/Pk does not diverge. Since pk does not exceed 
Pmax, we have the inequality 



Ail = V(l-e')/(A^-l), A^2=e/VN^, (15) 

for i 2 (the values of Aj^ for j > 2 do not enter into 
the analysis). Here e is a small number, and we will ig- 
nore terms of order e or higher. Now let pi = Pmax, and 
consider the measurement record where zero photons are 
detected in modes 3 to TV, and D photons are detected 
in output mode 2. To determine Cn^, note first that 
A22 ^ A2i for « ^ 2, so we may ignore those terms in the 
sum for Ss,n where A22 does not appear. Each term has 
magnitude A"]' A22A^j~^ 11], and there are D{D+ni-iy. 
such terms. Therefore, provided S2 = 1, 
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(16) 



s:Ss = D+l fc:5fc = l 



Here we are able to omit the condition P^ ^ because 
terms with Pg = are zero anyway. We may re-express 
the equation for ci as 
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We therefore obtain 



ci<K' E ^« E I'^'^i' 



s;S,=r>+l 



k;Sk=l 



Combining Eqs. (|ll|l and H13(l gives 
- < {M-D)^^ 

Co 1 - Pr 



(12) 
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This yields an upper limit on the ratio between the 
one and zero photon probabilities. One application of 
this result is that it is impossible to get one photon with 
unit probability, as it would be necessary for this ratio 
is infinite. Another consequence of (|14|l is that for D — 
M — 1 (i.e. the number of photons detected one less the 
maximum input number) an improvement can never be 
achieved. This case is important because it is the most 
straightforward way of eliminating the possibility of two 
or more photons in the output mode. 

In the remainder of this article, we will investigate sit- 
uations in which the single-photon contribution can be 
enhanced. As M < A^, and D > 0, the upper limit on 
the improvement in ci/cq is simply A'^. This is also the 
upper limit in how far ci can be increased above Pmax- 
We will now consider a scheme that gives a linear im- 
provement in ci/co, though not as high as N. In order 
to obtain a large value for the ratio ci/co, we want the 
inequality in Eq. H13|l to be as close to equality as possi- 
ble. In turn, this means that we want the vectors (An,) 



In the summation for c„j, we have different 
combinations of inputs such that Sg = D+ni and S2 — 1. 
Combining these results, we have 

"I „ l-Fmax i^max; 

^ (jV-l)!iJ^(i^ + ni-l)! ,2„,,2 ,213-2 

^"-11 ^'-22^'-21 
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(17) 



^iV-iy nil{N ~ D ~ ni)\' 

We have combined those factors that do not depend on 
ni into a new constant K", and used An w — 1. 

Using Eq. ifTTjl we find that 
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The maximum improvement in the ratio Ci/cq is obtained 
for D = \N/2], where Ci/cq = RlN^/A\/(N - 1). The 
multiplying factor [iV^/4j/(A^ — 1) is larger than 1 for 
all TV > 4. Thus we find that, provided there are at 
least 4 modes, we may obtain an improvement in the 
ratio ci/cQ. For p^ax < 1, ci « Pmax LA^V4J /(A^ - !)• 
For large N, the probability of a single photon increases 
approximately as 7V/4. This is linear with N, but is not 
quite as large as the upper limit of N. 

However, there are some drawbacks to this interferom- 
etry scheme [i.e. using the interferometer (|15() ]. The first 
drawback is the two-photon contribution. We find that 



C2/C1 {D + 1){N^D-1) 



ci/co 



2D(N - D) 



(19) 



which, for D — [iV/2], is less than 1/2, but it is close 
to 1/2 for large N . Since this is the same ratio as for a 
Poisson distribution, this is equivalent to using a coherent 
state. The two-photon contribution can be reduced by 
using larger D, but this is at the expense of reducing ci. 
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The second drawback is that improvements are only 
obtained for small Pmax- Although the improvement in 
the ratio Ci/cq is independent of Pmax, improvements in 
Ci can only be obtained for values of Pmax below 1/2. 
That is, this method can only be used to obtain im- 
provements in the probability of a single photon up to 
1/2, but not to make the probability of a single photon 
arbitrarily close to 1. Note that the above method only 
gives Ci > Pmax for four or more modes. We will now 
show that it is impossible to obtain an improvement in 
the probability of a single photon with fewer than four 
modes, and for various combinations of detections with 
larger numbers of modes. 

We first examine the case D = 0. Then we have only 
one term in the sum for cq, and cq = K'Pq. The expres- 
sion for ci becomes 

N N 

= K'RPo = coR. (20) 

Thus we have shown that ci/cq < R, so ci < Pmax- Hence 
there can be no improvement in the photon statistics if 
zero photons are detected. We can also obtain a similar 
result for the case D = 1, provided all the input pi are 
equal. In that case, we have 

co = K'J2 j^^^Po \A2kf = K'RPo. (21) 

K 

The value of Ci is given by 

k l;ljtk 

< ^K'R^Po J2 l^u^^k + AikA2if = K'R^Po. (22) 

k,l 

In the last line we have used the fact that Aik and 
are orthonormal. Thus we again find ci/cq < R, so ci < 

^^max* 

These results clearly eliminate the possibility of im- 
proving the probability of finding one photon with a two- 
mode interferometer. We have shown that detecting zero 
photons does not give an improvement, and if one photon 
is detected, then we must have M — £> = 1 or 0, so there 
again can be no improvement. Along the same lines we 
can also eliminate the three-mode interferometer. 

We have shown that it is impossible to improve the effi- 
ciency of a single-photon source by channeling more than 
one low-efficiency single-photon state into a linear optic 
interferometer and detecting all but one of the photons. 
This eliminates the most straightforward scheme for ob- 
taining an output state with no more than one photon. 



It is possible to obtain an improvement for more general 
detection results, but at the expense of non-zero prob- 
abilities for two or more photons. We have not proven 
that it is impossible to obtain an improvement in the 
probability of a single photon while restricting to zero 
probability for two or more photons; however, numerical 
searches indicate that it is unlikely. 
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